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Abstract. The purpose of these notes is to give a brief review of superfluidity in 
neutron stars. After a short presentation explaining why and how superfluidity is ex- 
pected in the crust and core of neutron stars, consequences on thermal evolution and 
rotational dynamics are discussed. The second part summarizes a formalism that has 
been recently developped to describe the hydrodynamics of superfluids or supercon- 
ductors in the framework of general relativity. As an application, one can compute the 
oscillations of a two-component relativistic neutron star. 



1 Introduction 

The present contribution will discuss superfluidity and superconductivity in a 
rather extreme environment: neutron stars. Neutron stars, dense stars composed 
mainly of neutrons, were envisaged by Landau as soon as the neutron was dis- 
covered in 1932. In 1934, Baade and Zwicky suggested that supernovae were 
manifestations of a transition from an ordinary star to a very dense neutron 
star. In 1967, Hewish and Bell discovered the first "pulsar" (pulsating source of 
radio), which was soon identified as a rotating neutron star. The following year, 
two new pulsars, much studied until now, were discovered: the Vela pulsar with 
a period P = 89 ms and the Crab pulsar with P = 33 ms. Since then, more than 
one thousand pulsars have been detected. 

Neutron stars are rather impressive objects. They contain a mass of the order 
of the solar mass confined in a radius i? ~ 10 km, which implies an average mass 
density of the order p ^ lO^'^g/cm'^ . They can rotate up to several hundred 
times per second. Due to these extreme conditions, neutron stars are of interest 
for various branches of physics. First, neutron stars are so dense and so compact 
that their gravitational field is very strong, 

and can be described correctly only with general gelativity. Moreover, binary 
pulsars have been extraordinary "laboratories" to test the strong field predic- 
tions of general relativity, in particular to verify (indirectly) the existence of 
gravitational radiation. 

Second, the magnetic field, typically of the order of 10^^ G, plays a very 
important role in the neutron star physics, in particular in the pulsar emission 
mechanism (see e.g. |1[). Third, neutron stars are unique places where one can 
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find matter in an ultra-dense state. The density is indeed higher than the atomic 
nuclear density. The corresponding pressure therefore depends essentially on 
the strong interactions, and neutron star observations are potentially a very 
rich source of information about the behaviour of strong interactions at high 
densities, which, until now, remains very poorly known from both theoretical 
and experimental points of view. 

Finally, neutron star physics has also some connection with the low tempera- 
ture physics. Studying the nucleon interactions, Migdal suggested the possibility 
that neutron star matter become superfluid at sufhciently low temperature. If 
the temperature of neutron stars, typically T ~ \Q^K, may seem huge in ter- 
restrial standards, it is in fact smaller than the superfluid critical temperature 
evaluated as ~ 10^ - 10^° K. 

In these notes, we will focus on two aspects of the physics of neutron stars: 
the property of superfluidity with its consequences on the evolution of a neutron 
star; the necessity to consider a neutron star as a general relativistic object. 
Many other interesting aspects of neutron stars will thus be left aside and the 
curious reader will find some useful information in e.g. 

The plan will be the following: we will begin with some details about the 
superfluid properties of the interior of neutron stars. Observable consequences 
of superfluidity will then be discussed: first, the impact on the thermal evolution 
of neutron stars, then on their rotational dynamics. In a second part, we will 
summarize some elements of a formalism that gives a hydrodynamical description 
of the interior of neutron stars in the framework of general relativity. Finally, 
as an illustration, oscillations of two-component superfluid neutron stars are 
condidered. 

2 Superfluidity and superconductivity 
2.1 Composition of the interior of a neutron star 

Let us start by describing the matter composition inside a neutron star. Because 
there is a strong density gradient from the exterior to the interior of the star, 
the composition changes dramatically with the radial distance from the center. 
Let us progress from the surface of the star towards its interior. 

Ignoring here the surface ocean, the first layer of a neutron star is the outer 
crust: it is made of a lattice of nuclei with a gas of relativistic degenerate elec- 
trons. The nuclei are richer and richer in neutrons as the density increases (see 
e.g. 1^ for a recent review). At some point {p ~ 4.3 x lO^^g/cm'^), neutrons 
begin to leak out of the nuclei. This is called the neutron drip transition and this 
marks the boundary between the outer crust and the inner crust. For the latter, 
one still has a lattice of nuclei immersed in an electron gas, but in addition there 
is a liquid of free neutrons. The density still increasing, the difference between 
the neutrons in nuclei and the neutrons outside becomes fainter and fainter until 
the nuclei simply dissolve (at a density p ^ 2 x 10^'^g/cm^). One has reached 
the limit of the crust and one now enters the realm of the core. 
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In the outer core, one will find the coexistence of a gas free neutrons, which 
is by far the main species, with a liquid of free protons, plus a gas of electrons 
(and muons) which ensures charge neutrality. Finally the deeper layers of a 
neutron star, the inner core, remain mysterious. Several possibilities can be 
envisaged, including appearance of hyperons, condensation of pions or kaons, or 
the transition from hadronic to quark matter. 

2.2 Energy gaps and critical temperature 

Like electrons in a superconductor, the superfluidity of neutrons in neutron stars 
is due to the pairing of two neutrons near the Fermi surface in momentum space, 
according to the Cooper mechanism (see e.g. 0). In the inner crust, the neutron 
Cooper pairs are preferably in a state ^5*0. In the outer core, where the density 
is much higher, the neutron Cooper pairs are of the '^P2 type (see § and Q). 
The evaluation of the energy gaps give typically values of the order of 1 MeV. 

In the same layer, the protons are free and can undergo the same process 
as neutrons, i.e. form Cooper pairs. Being fewer than the neutrons, they are 
expected to condense into a ^5*0 state 0. Note that hyperons that may exist in 
the inner core are also expected to be superfluid. 




inner crust 



outer crust 



Fig. 1. Schematic view of the internal structure of a neutron star 
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2.3 Various equations of state 

In order to determine the global structure of neutron stars, one needs the effective 
bulk equation of state of neutron star matter. Inserting this equation of state 
in Einstein's equations, it is immediate, in the case of a non-rotating star (see 
subsection 6.1), to obtain the radial profile of the star, i.e. the radial profile of 
the energy density, of the pressure and of the metric coefficients. One can also 
take into account the rotation of the star. 

The main problem however is that the equation of state is unknown at very 
high densities. The reason is the lack of experimental data as well as the theo- 
retical and computational challenge to evaluate the interactions of high density 
matter. Therefore there exist many different equations of state in the literature 
(see e.g. ^ for a recent review) and the hope is to infer from observations of 
neutron stars some constraints on the high density equations of state (see Q 
and 0). 



3 Cooling processes in neutron stars 

Born with a temperature of the order 10^^ — 10^^ K, neutron stars cool very 
rapidly to temperatures of less than 10^° K within minutes. The subsequent ther- 
mal evolution is more uncertain and could be strongly affected by the presence 
of superfluidity. The mechanism by which neutron stars cool down is essentially 
neutrino emission. Neutrinos can be produced in two types of processes: 

• Direct URCA processes 

They correspond to the simplest beta reactions 

n ^ p + e~ + De, p + ^ n + Ve- (2) 

• Modified URCA processes 

They correspond to beta reactions with the presence of a second nucleon. 

Direct URCA ractions are (not surprisingly) more efficient than the modi- 
fied ones but they can occur only if conservation of energy and momentum is 
satisfied, which implies some constraint on the species fractions, namely that 
the proton fraction is sufficiently high. If this is not the case, then direct URCA 
processes are completely suppressed and cooling occurs only through modified 
URCA processes and is much slower. 

However, taking into account the existence of superfluidity leads to the con- 
clusion that the direct URCA processes could be partially suppressed: the neu- 
trino emissivity would be reduced by a factor exp(— Z\/fcT), where A is the 
superfluid gap energy. The reason is that, before the beta reaction takes place, 
one needs to break the Cooper pair. 

The existence of superfluidity therefore enables us to envisage, between the 
two extreme scenarios suggested above of fast cooling or slow cooling, an inter- 
mediate scenario of moderate cooling, which may be testable in the future by the 
temperature measurements of neutron stars (see [ll| ] and references therein). 
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4 Rotational dynamics of neutron stars: glitches 

Another important consequence of superfluidity concerns the rotational dynam- 
ics of neutron stars, which can be followed extremely precisely via the observation 
of the radio signal from pulsars. The main feature is an extremely stable peri- 
odic signal, the pulsar being analogous to a lighthouse beacon: the radio emission 
is indeed coUimated along the magnetic axis which does not coincide with the 
rotation axis. 

For isolated pulsars, one observes a steady, although tiny, increase of the 
pulse period, i.e. a steady decrease of the angular velocity. This is interpreted as 
a loss of angular momentum due to the emission of electromagnetic radiation. 
Assuming a magnetic dipole, the rotational energy loss is given by 

E^wn = — , 3 

where I is the moment of inertia and 9 the angle between the magnetic and 
rotation axes. 

In addition to this slow increase of the period, a few pulsars exhibit some 
rare events named glitches at which the period suddenly decreases, i.e. the 
angular velocity suddenly increases, followed by a slow exponential relaxation 
with a typical time scale of the order of days to months. The two most famous 
'glitching' pulsars are the Vela and Crab pulsars with period changes Af2/f2 of 
the order of 10~^ and 10~* respectively. 

The current explanation for these glitches is the following. It is based on the 
idea that the neutron star interior contains a weakly coupled component that 
is not directly slowed down by the electromagnetic torque and therefore rotates 
slightly faster than the star crust. What is directly observable is the angular 
velocity of the crust, and if, for some reason, there is a rapid transfer of angular 
momentum from the faster component to the crust, then one should see a sudden 
increase of the pulsar period. Of course, a natural candidate for the component 
rotating faster is the neutron superfluid. 



4.1 The tv^ro-fluid model 

The relaxation after the glitch can also be explained, at least qualitatively, by the 
simple two-component model of Baym et al . One considers two components 
with respective angular velocities f2c and and respective moments of inertia, 
Ic and /„. The first component corresponds to the crust and to whatever is 
strongly coupled to it, whereas the second component corresponds to the neutron 
superfluid weakly coupled to the crust. Note that the pulsar angular velocity I? 
which is the directly observable quantity can be identified with the crust angular 
velocity, so that f2c ~ ^- The evolution of the two components is governed by 
the system 



/cl? = — -{SI — fin) — oi, 

Tc 



(4) 
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Inf2r, = -{0-nn), (5) 



where the first term on the right hand side represents a coupling between the 
two components, couphng characterized by the timescale Tc- a corresponds to the 
electromagnetic torque. One can integrate the above system just after a glitch 
characterized by a sudden jump of the crust angular velocity Af2o = fl{t = 
0+) — n{t = 0~). The subsequent evolution is then 

n{t) = f2o (t) + AHo (Qe-*'^ + 1 - q) , (6) 

where Q is the healing parameter and f2Q{t) = J?o — {a/I)t (J?o being a constant) 
corresponds to the evolution without glitch, and the characteristic relaxation 
timescale is given by 

r = |re, (7) 

with 7 = 7„ + 7c. 



4.2 Role of the vortices 

The two-fluid model gives an effective view of the rotational evolution, but it is 
interesting to explore in more details the underlying mechanism responsible for 
the coupling between the two components. This is where the vortices enter into 
play. 



1 l i 




t 

Fig. 2. Evolution of the angular velocity during a glitch 
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The neutron star at its birth is a rapidly rotating object. With its coohng, the 
temperature reaches the critical temperature under which superfluidity appears. 
However, one essential property of a superfluid is that its flow is irrotational. The 
way the superfluid solves the contradiction between irrotationality and rotation 
of the star is, like superfluid Helium in laboratories, via the creation of quantized 
vortices that carry all the angular momentum of the superfluid. The superfluid 
'velocity' (in fact, momentum) being the gradient of a phase, 

Pn = 2mnVn = ?iV<^, (8) 

the circulation around any closed path is given by 

K = \ dlvn = {2ttN), (9) 

where N is an integer. Superfluid vortices usually correspond to one unit (N = 1) 
of quantized circulation n = h/ (2m„). One can then easily compute the average 
density of vortices for a superfluid rotating at uniform angular velocity f2„: 

ny = 4^^. (10) 



Putting numbers, one finds 



nv = 6.3 X 10'^ \~) '^o'^^i*^^^ P^"^ (-'■-'-) 

In the neutron star core, the existence of a proton superconductor, which 
is believed to be a type II superconductor, suggests the presence of magnetic 
vortices, or fluxoids, which carry the magnetic flux going through the neutron 
star. Note that one must be careful to compute the total energy of a fluxoid in 
a rotating superconducting background JTsf . Moreover, because of the 'entrain- 
ment effect' by which the neutron superfluid momentum is associated with both 
neutron and proton currents, neutron superfluid vortices, in the neutron star 
core, will carry a magnetic flux as well, thus implying a strong coupling between 
the core superfluid and the crust pUTsI] . 

Now, the vortices play an essential role for the coupling of the superfluid to 
the normal component of the star, because the pure superfluid has no viscosity 
at all and, as such, can circulate without affecting any normal component. But 
the vortices interact with the normal component. Without entering the details, 
one can distinguish two types of interaction between the vortices and the normal 
part of a neutron star: 

• Pinning 

This is the situation where a vortex, or a segment of it, is anchored to the 
(inner) crust. The pinning is due to the interaction between the vortex and 
the nucleus, which can be attractive or repulsive depending on the matter 
density 0,0 
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• Friction 

This is the case where the vortex, moving through the normal component, 
interacts with it in some dissipative processes, which results into an effec- 
tive friction force (per unit length) exerted on the vortex by the normal 
component of the form 

Fd,^,g^C{vv -V,). (12) 

Let us assume here that we are in the friction case. The vortices, interacting 
with the normal component, are also sensitive to the motion of the superfluid 
through a Magnus force term. 



' M 



iK X {W - Vn), (13) 



which is a force per unit vortex length. The effective motion of the vortices is 
obtained by requiring that the two forces (|l^) and ( [l^ ) just cancel. This implies 
that the vortices must have a radial motion in addition to the angular motion. 
One finds easily that the angular motion is given by 

Hv = ''r'^n+''-^- ^ (14) 
Cr ~\~ Cr 

which means that the angular velocity of the vortices is simply a weighted average 
of the superfluid and normal angular velocities, the weight depending on the 
dimensionless friction coefficient 

Cr = — . (15) 

The two extreme cases are: the case where the friction is zero, the vortices being 
then in corotation with the superfluid; the case where the friction is huge, the 
vortices being then in corotation with the normal component. As for the radial 
velocity of the vortices, it is given by 



Cr ^ Ci 



-r. 



(16) 



One can then compute the effective coupling timescale between the two com- 
ponents of the stars, by noting that the evolution of the angular velocity of the 
superfluid component is directly given by the radial velocity of the vortex array, 
according to the expression 

= -2^/2„, (17) 
r 

simply because the angular momentum of the superfluid is proportional to the 
vortex density. Let us consider now the evolution of the angular velocity lag 
between the two components, 

w = i7„ - I2„ (18) 



Superfluidity in relativistic neutron stars 9 



and let us define the coupling timescale as = Combining the expressions 

obtained above, one finds 

The above Newtonian analysis has been generalized to a general relativistic 
context in iQ. The coupling between the crust and the superfluid can also be 
evaluated when vortices are pinned, in which case vortices can move outwards 
by vortex creep, the temperature being a crucial parameter [ p^ . 

4.3 Origin of the glitches 

There is no certainty at present on the physical origin of observed glitches. The 
scenario which seems to have attracted the greatest attention is the one proposed 
by Anderson and Itoh ||20|] . In this model, vortices are supposed to be pinned in 
the crust due to the nucleus- vortex interaction. Because the normal component 
is slowing down, a Magnus force, due to the superfluid motion relative to the 
vortices, progressively builds up until this force reaches a critical amplitude for 
which the pinning sites simply break. The vortices then can suddenly move 
outwards, thus transfering angular momentum from the superfluid to the crust, 
which explains the sudden spin- up of the crust. 

Other models have been proposed to explain glitches. For example, a glitch 
could be induced by a sudden heat wave propagating in the star, which would 
increase the effective coupling between the superfluid and the crust |2^. Such 
a heating could be generated for instance by a crustquake. In this respect, one 
must mention that the differential rotation between the superfluid and the crust 
generates a centrifugal buoyancy force which might increase the crust stresses 
and maybe lead to a glitch [ p2[ . 

5 Relativistic description 

Until very recently, only Newtonian theory was used in the studies of superflu- 
idity in neutron stars, whereas, in parallel, general relativistic studies of neutron 
stars, even numerical, were based on perfect fluid matter. The purpose of the 
work summarized here is to establish a bridge between these two approaches, 
taking into account both the highly relativistic nature of neutron stars and their 
superfluid interior. We shall present a macroscopic formalism ||2^] allowing for 
the average effect of vorticity quantisation in a rotating superfluid, which can 
then be extended j2j] to describe as well superconducting fluids, such as pro- 
tons in the core of neutron stars. The Newtonian version of this formalism can 
be found in [^ . In the next section, a simplified relativistic two-component 
model of neutron star, which allows for differential rotation of the superfluid 
component, will be used to study oscillations of neutron stars. 
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5.1 Perfect fluid in general relativity 

Let us start by recalling the simplest case, that of the perfect fluid. In general 
relativity, a perfect fluid is characterized by 

• a four-velocity vector (where /i = 0, 1, 2, 3 is a spacetime index denoting 
the time coordinate and the three space coordinates), which satisfies the 
normalization gp^u^'u'^ = — 1 where gp^y is the metric tensor describing the 
geometry of spacetime (and which is used to raise or lower the indices, e.g. 
Up = gpau"). 

• a particle number density scalar field n{x^) giving in each spacetime point 
the density of particles as seen by an observer comoving with the fluid. 

• an equation of state of the form p = p{n) giving the energy density p as a 
function of the particle number density n. 

Note that we have chosen here the simpler case of a barotropic perfect fluid 
where the equation of state depends on only one parameter. Note also that one 
could characterize the fluid by its energy density instead of its particle number 
density, and then derive the number density by inverting the equation of state. 

The equations of motion for the relativistic fluid will be generalizations of 
the well-known Newtonian equations of fluid mechanics, namely, 

• the matter conservation equation, which can be written in the present con- 
text in the very simple form, 

V^n'' = 0, (20) 

where n'^ = nu'^ is the particle current and stands for the covariant 
derivative. 

• the relativistic Euler equation, which can be written in the very condensed 
form 

n''^[p^la] = 0, (21) 

where the brackets represent antisymmctrization on the indices and where 
one has conveniently introduced the momentum covector p,p = pUp, p, being 
the relativistic chemical potential deflned by 

It is possible to derive directly the above equations of motion from a variational 
principle, the Lagrangian density being simply the energy density considered as 
a function of the particle number current n^. Any variation of is not allowed, 
but only the convective variations, which correspond to variations of the particle 
flow lines (see for more details). 
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5.2 Relativistic superfluid 

What will characterize a "pure superfluid" (i.e. a superfluid at zero temperature, 
without gas of excitations), with respect to the more general class of perfect flu- 
ids, is the fact that its motion is locally irrotational. In the standard language, 
this is a consequence of the fact that the 'superfluid velocity' can be expressed 
as the gradient of a quantum scalar phase. Strictly speaking, the standard 'su- 
perfluid velocity' is in fact a momentum divided by a somewhat arbitrary mass. 
In the relativistic context, this will be generalized by the expression 

/ia- = ?iVo-(/?, (23) 

so that the property of irrotational flow can be written as the vanishing of the 
vorticity tensor 

Wp,=2Vyp^i„y (24) 

The previous condition is valid locally. However, as explained above, a su- 
perfluid can exhibit a non-irrotational flow at the price of being threaded by 
quantized vortices. This means that macroscopically, on distances bigger than 
the typical intervortex separation, one would like to describe the superfluid as 
well as the vortices in an average way. This can be done by considering the 
vorticity tensor Wpo-, now non-zero since there are vorticies, as a fundamental 
quantity. In fact, this tensor contains information about the density of vortices 
as well as about the local average direction of the vortex array. 

The idea then is to construct a generalized Lagrangian density that depends 
not only on the particle current n'^ but also on the vorticity tensor Wpa. The 
variations with respect to these fundamental variables will define the canonical 
momenta, 

5A = fi^Sn" - ^VSwp^. (25) 

This can be seen as a generalization of the perfect fluid where A ~ —p and only 
the first term is present on the right hand side. 

The equations of motion are the matter conservation equation (|2^) , as in the 
perfect fluid case, and a generalized Euler-type equation of the form, 

WpAn" - VpA"") = 0. (26) 

Of course, the expression for A'"^, which essentially represents the energy density 
and momentum of an individual vortex, must be specified by resorting to a 
'microphysical' model of the vortices (see |^ and |2|] for relativistic descriptions 
of vortices) . 

5.3 Superfluid-Superconducting mixtures 

It is useful to extend the previous formalism to the case where the superfluid 
particles (or Cooper pairs) are electrically charged, like electrons in laboratory 
superconductors, in order to be able to describe the protons in the cores of 
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neutron stars. Let us consider several species labelled by the index X, with 
respective particle currents and respective electric charge per particle e^. 
The total electric current is then given by 

/ = E^'<- (27) 

X 

Proceeding as before, In order to obtain the equations of motion for the 
global system of fluids, one uses a variational principle based on a Lagrangian 
density C that will be the sum of three contributions: 

• a "matter" contribution tIm, which depends only on the 'hydrodynamical' 
part of the system, and which is a function only of the particle currents 
nP, or more exactly of all scalar combinations obtained by their mutual 
contractions, 

• an electromagnetic interaction term depending on the electromagnetic gauge 
form Ap, 

Ai = fAp , (28) 

• an electromagnetic field contribution 

= -;^Fp^F^P . (29) 

where the electromagnetic field tensor Fp^- is defined in terms of the gauge 
form by the usual expression Fp^ = 2^[pA^]. 

Considering the variations of £ = Am + Ay + Ai with respect to the particle 
currents naturally suggests to define the canonical momentum covectors 

TT-p = fi-p + e-Ap (30) 

which are the sum of a pure 'hydrodynamical' (or 'chemical') part and of the 
electromagnetic gauge form weighted by the electric charge of the species. From 
this, it is convenient to define generalized vorticity tensors defined by 

w-p, = 2V[p^^] = 2V[,M\] + e'^Fp^ . (31) 

The equations of motion for the system consist of the separate matter con- 
servation equation for each species (one can generalize to allow for chemical 
reactions between various species |l^), 

VpnP = 0, (32) 

and of Euler-type equations, which can be written in the very compact form 

n^w-p = . (33) 

Note that in the case of a charged component, this equation contains the electro- 
magnetic field tensor and thus automatically includes for instance the Lorentz 
force exerted on the fluid. 
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It is also useful to write for the system under investigation the associated 
total stress-energy-momentum tensor, which appears on the right hand side of 
the Einstein's equations, 

- = 8ttGT^,, (34) 

if one needs to determine the spacetime metric. Once again, the variational 
principle is very useful because variation of the Lagrangian density with respect 
to the metric directly yields 

where the part derived from the material Lagrangian density contribution tIm 
is given by 

X X 

while the electromagnetic contribution has the usual Maxwellian form 

V. = ^ (P'^'F'^'^ ~ {p^-F.^gP^) . (37) 

The above relations are valid for ordinary perfectly conducting fluids. The 
condition characterizing superconductors, generalizing the superfluid condition 
(p3|), is that the generalized momentum should be the gradient of a quantum 
phase, 

TT- =M-, + e- =Wp^-, (38) 

so that superconductors (including the case of superfluids) are characterized 
locally by a vanishing generalized vorticity tensor. Of course, it is then possible, 
in analogy with the treatment of the previous subsection, to extend the formalism 
in order to describe on average a superfluid-superconducting mixture threaded 
by superfluid vortices as well as magnetic flux tubes. 

6 Relativistic neutron stars 

Let us now see how one can construct a concrete model of neutron star in a 
general relativistic framework. To simplify, we will assume that the neutron star 
is made of simply two components, each being treated as a perfect fluid. First, 
we will consider the case of a static, i.e. nonrotating, neutron star. Then we will 
consider small, i.e. linear, oscillations of the two-component neutron star. 

Let us call, loosely, the component corresponding to the free neutron super- 
fluid the "neutron" component and the component corresponding to the normal 
fluid the "proton" component. Their respective particle currents will be denoted 

= nuP and pP = pv^ , where and are unit four- velocities, which are 
not aligned if the two components are not comoving. Their equations of motion 



consist of two separate conservation equations of the type (20) and of two Eulcr 
equations of the type (|2l|). Here, the two components are treated as ordinary 
perfect fluids and one ignores the refinements due to the presence of vortices in 
the superfluid component. 
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6.1 Static neutron star 

We first consider the configuration corresponding to a non-rotating neutron star, 
which is thus static and spherically symmetric and for which the metric can be 
written in the form 

ds^ = g^^dz'^dx" = -e'^^^'Ut^ + e^^^'Ur^ + [dO^ + sin^ ed<t)'^) . (39) 

The particle currents, since the two fluids are motionless, are necessarily of the 
form 

n^ = nM^, p^'=pu^', = {u", 0, 0, 0}. (40) 

Writing Einstein's equations with the above ansatze, one obtains the well-known 
TOV (Tolman-Oppenhcimer-Volkov) equations: 

dP _ {p{r) + Pir)){mir)+Anr^P{r)) K^M^p^.d!. 
H^--^ r2(l-2Gm(r)/r) " " 2^^^"^^ + ^^'^^ d^^' 



and 



eAM=(l_2^') , (42) 



with m{r) = AttG dr' p{r')r'^ . The mass of the star is given by M = m(r), 
where R is the star radius defined by the condition that the pressure vanishes. 



6.2 Oscillations of superfluid neutron stars 

Oscillations of relativistic stars have been studied by many authors in the case 
of a single perfect fluid. Here, we summarize the first investigation for a two- 
component relativistic star p9|, thus generalizing two previous studies of two- 
component Newtonian stars ]30| . 

In a relativistic approach, one must consider not only the perturbations of 
the quantities describing the matter, for example the particle density and the 
velocity, but also the perturbations of the spacetime metric 17^^. After decom- 
position of the perturbations into spherical harmonics, labelled by I and m, it 
is convenient to distinguish even-parity (or polar) perturbations and odd-parity 
(or axial) modes depending on their transformation under parity. We will re- 
strict ourselves here to the case of even-parity perturbations, which is the most 
interesting since scalar quantities, such as the matter densities, give only this 
type of perturbations. Also, one needs to study only m = perturbations, since 
there is a degeneracy in m for each given I because of the spherical symmetry of 
the background. 

Choosing a specific gauge, one will consider metric perturbations of the form 

Sgoo - -e-'r'Hoe'^'Piie), 6gor = 5gro = -iiur'+^ H^e'^' Pi{e), 

Sgrr - -e^r'H^e'^'Piie), 5gee = Sg^^./ sin^ 6 = -r'+^ Ke'^' Pi{9), (43) 



Superfluidity in relativistic neutron stars 15 



where Pi stands for the Legendre polynomial of order /. The matter perturbations 
will be described, in a Lagrangian way, by the matter displacements, 

Cp = r'-^e-^/^Wpe^^'Piie), = -r'-^VpC^'^'dePiie), (44) 

from which one can compute the perturbed velocities, by taking the time deriva- 
tive or the variations Sn and Sp by using the perturbed conservation equations. 
Inserting the above expressions for the perturbations into the perturbed Ein- 
stein's equations as well as the perturbed Euler equations, one ends up with a 
system of linear equations, consisting of two contraints (one is H2 — Hq, the 
other expresses Hq in terms of the other perturbations) and of first order differ- 
ential equations of the form 

^ = Qi.^Y, (45) 
ar 

where Y is 6-dimensional column vector containing Hi, K, Wn, Vn, Wp and Vp, 
and Qi_cj is a 6 x 6 matrix with r-dependent coefficients that depend only on 
the background configuration, as well as on / and uj. Considering the boundary 
conditions at the center and at the surface of the star, this system can be solved, 
up to a global amplitude, for any value of lo. The physically relevant modes, 
however, also called quasi-normal modes, correspond the specific values of lo for 
which the metric outside the star represents only outgoing gravitational waves. 

A numerical investigation for a very crude, and unrealistic, model of two 
independent polytropes (with different adiabatic indices) has shown that the 
two-component star will exhibit new modes, superfluid modes, which are specific 
of the existence of two components since the two fluids are counter-moving for 
these modes. 
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